Landau- Zener Interference in Multilevel Superconducting Flux Qubits Driven by 

Large Amplitude Fields 
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We proposed an analytical model to analyze the Landau-Zener interference in a multilevel su- 
perconducting flux qubit driven by large-amplitude external fields. Our analytical results agree 
remarkably with those of the experiment [Nature 455, 51 (2008)]. Moreover, we studied the effect 
of driving-frequency and dephasing rate on the interference. The dephasing generally destroys the 
interference while increasing frequency rebuilds the interference at large dephasing rate. At certain 
driving frequency and dephasing rate, the interference shows some anomalous features as observed 
in recent experiments. 

PACS numbers: 74.50.+r, 85.25.Cp 



I. INTRODUCTION 

Superconducting Joscphson devices coherently driven 
by external fields provide new insights into fundamen- 
tals of quantum mechanics and hold promise for use in 
quantum computation as qubitsi 1 -' 2 -' 3 - Compared with nat- 
ural atoms and molecules, these devices can be strongly 
coupled to external radio-frequency (rf) fields while pre- 
serving quantum coherence^ 5 -! 6 - A large number of exper- 
iments associated with strong driving have been done on 
these devices such as Rabi oscillations in the multipho- 
ton regime,—^ Bloch oscillations,— Landau-Zener(LZ) 
interferenc o 11 ' 12 i 13 and dressed states of superconducting 
qubits under extreme driving.— 

Recently, coherent dynamics of superconducting qubits 
in the regime dominated by LZ transitions were exten- 
sively studied, 11 ' 12 - 13 ! 15 ^ 6 ! 17 - 18 ! 19 ! 20 ! 21 ! 22 In this case, the 
driving frequency is much smaller than energy-level sep- 
aration and the transitions occur via the LZ process at a 
level crossing . 15 ' 16 One may use LZ transitions to enhance 
the quantum tunneling rate ; 18 ' 19 to prepare the quantum 
stated to control the qubit gate operations effectively^, 
and so on. Moreover, repetition of the LZ transitions 
can induce quantum mechanical interference, which leads 
to Stueckelberg or Ramsey-type oscillations ! 23 ' 24 Since 
the theoretical scheme of observing LZ interference in 
qubits was proposed by Shytov et alji 5 . a series of beauti- 
ful experiments on LZ interference were implemented in 
two level systems such as flux qubit a 11 ' 13 ' 17 ' 22 and charge 
qubits, 12 which provided an alternative method to ma- 
nipulate and characterize the qubit in the strongly driven 
regime. 

A new regime of strong driving was reported in a recent 
worki^ 5 . Unlike previous experiments which employed a 
two level system ; 11 ' 12 ' 13 the qubit in this experiment was 
driven through a manifold of several states spanning a 
wide energy range. The population of the qubit under 
large-amplitude fields exhibited a series of diamond-like 
interference patterns in the space parameterized by flux 
detuning and microwave amplitude. The interference 
patterns, which displayed a multiscale character, encoded 



the information of several energy levels of the system. In 
a recent work by Rudner et al.^ 6 - a skillful and concise 
method employing a two-dimensional Fourier transform 
was used to study such a system. They transfer the mea- 
surement results of the energy domain to the time do- 
main, and realize a "tomogram" of the time evolution of 
the qubit phase, from which the decoherence time of the 
qubit can be easily obtained. 

In this article, rather than extract information by de- 
composing or translating the interference patterns, we 
start from the point of reconstructing the interference 
patterns through analyzing the dynamics of the system, 
i.e., we map the time variable (including both decoher- 
ence rate and field-driving rate) onto the distribution of 
the qubit 's population. Our model used to analyze mul- 
tilevel systems are based on the well developed theory of 
two-level systems ! 13 ' 15 ' 27 ' 28 ' 29 ' 30 This article is organized 
as follows. In Sec. II the basic models were introduced 
firstly, from which we used rate equations to discuss the 
dynamics of the two diamond-like interference patterns 
separately. The results agree with those of experiments 25 
very well. In Sec. Ill, we discussed the effect of driving 
frequency and dephasing rate on the LZ interference. By 
tuning the driving frequency and dephasing rate one can 
expect different interference patterns. Some anomalous 
interference patterns, such as the moire-like pattern re- 
ported in a recent experiment ^ can be well explained 
using our model. 



II. MODEL AND METHOD 

We start from a driven two-level system subject to the 
effects of decoherence; 13 ' 27 



H(t) 



A „ h(t) A 



(1) 



where A is the tunnel splitting; a x and a z are Pauli ma- 
trices. h(t) is the time dependent energy detuning from 
an avoided crossing: 



h(t) = e + Asintot + 8e nmse (t), 



(2) 
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where e is the dc component of the energy detuning, A 
and li are the amplitude and frequency of the driving 
rf field respectively, Se no i se (t) is the classical noise. As 
discussed in ref [13], by using white noise model and per- 
turbation theory, the rate of LZ transitions between the 
states |0) and |1) can be obtained: 



W(e,A) 



A 2 



E 



n 



(3) 



where T2 = I/T2 is the dephasing rate and J n (x) are 
Bessel functions of the first kind with the argument 
x = A I 'li. Eq.Q implies that the transition rate is pro- 
portional to A 2 which is decided by the energy structure 
of the system. 

Extending Eq.Q to multilevel systems, the LZ transi- 
tion rate between states i) and \j) can be written as: 



Wi, 



i>A) 



E 



(4) 



where Ay is the avoided crossing between states |z) and 
\j), and eij is the dc energy detuning from the corre- 
sponding avoided crossing Ay . Eq. (0} can be derived in 
the same way as in ref [13] considering only the direct 
coupling between states \i) and \j). We emphasize that 
the situation will change for no decohcrcncc, in which 
coherence evolution among all coupled levels should be 
considered. While in our case, the coherence evolution 
is reduced to rate equations as discussed below, which is 
appropriate in dealing with stationary population distri- 
bution. 

Hereafter, we focus on the multilevel superconducting 
flux qubit, a superconducting loop interrupted by three 
Josephson junctions. Near flux bias $ = 0.5m$o, the 
system exhibits a double-well potential parameterized by 
the flux detuning (5$ = $ — 0.5m$o- Fig. 1 shows the 
lowest four levels of the system as a function of flux de- 
tuning. Whenever the diabatic states \i) (i —0, 1; right 
well with negative slope ) and \ j) (j =2, 3; left well with 
positive slope) are degenerate, avoided crossing Ay forms 
because of the inter-well tunneling. If a microwave flux is 
applied (Fig. 1, blue sinusoid), the detuning flux of the 
system is 

<5$(t) = 6$ dc + S^ac = 5<S>dc + Qrfsinut, (5) 

with $ r /, in the unit of flux, proportional to the mi- 
crowave source voltage. Then the time dependent energy 
detuning (see Eq.®) of states \i) and \ j) of the flux qubit 
can be expressed as: 



iy(t) = (|m; 



\mj\) {5$ dc + ®r f si-nut) + 8e noise (t), (6) 



where m, = dEi(&)/d<& is the diabatic energy-level 
sloped of state \i) in units of frequency per flux. From 
Eq.(d]) and Eq.©, the LZ transition rate Wjy is a func- 
tion of 5$rf c and $ r /, which can be easily controlled in 
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FIG. 1: (color online) Energy- level diagram of a supercon- 
ducting flux qubit illustrating the level-crossing positions for a 
particular static flux detuning 8§d c . Blue solid line represents 
the time-periodic detuning & r fsinujt. The red path and the 
green path represent two different transition processes leading 
to different interference patterns. The dotted lines represent 
intra- well relaxation processes. 



experiments. To understand simply, the tunnel split- 
ting Aij serves as a channel that connects different en- 
ergy states i) and \j) when the systems are driven by 
large-amplitude fields. Whether the channel is 'open ' 
or 'closed' is related with both the energy detuning ey 
and the amplitude of the driving field A. As observed in 
experiments, the qubit population exhibited a series of 
'diamonds' patterns 2 ^ in the space parameterized by $ r / 
and 6$>dc ■ Here we mainly discuss the first two diamonds 
produced by the transitions among the lowest four energy 
level states. The same methods can be used to analyze 
other 'diamonds' resulting from transitions among higher 
energy levels. In order to have a clear physical picture 
of the interference patterns, we analyze the first two dia- 
monds in part A and part B separately, and make a more 
quantitative analysis in part C. 



A. THE FIRST DIAMOND 

As shown in Figure 1, the red path represents the 
main transition process generating the first 'diamond'. 
In order to describe the time evolution of the popula- 
tion in the presence of driving fields, we employ a rate 
equation approach, in which the qubit level occupations 
Pi(i = 0,1,2,3) obey: 

Po = -PoWq2 +PiT w +P2(Wq2 +r 20 ), 

pi = -pi(Wi2 + r 10 ) +P2W12, 

Po+Pi+P2 = l, (7) 
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FIG. 2: (color online) (a) LZ transition rate W02 versus flux 
detuning and driving amplitude. (b)LZ transition rate W12 
versus flux detuning and driving amplitude, (c) Qubit pop- 
ulation in the left well obtained from Eq.©. The left edge 
and the right edge mark the parameter values where the level 
crossings A02 and A12 are first reached respectively. The 
parameters we used are from experiments 25 with driving fre- 
quency lu/2tt = 0.16GHz, A 02 = 0.013GHz,Ai2 = 0.09GHz. 
The locations of the level crossings A02 and A12 are Om"l>o 
and 8.4m<t>o respectively. The diabatic energy-level slope 
|mo|(|m2|)=1.44GHz/m$ , and |mi|(|m 3 |) =1.09GHz/m$ . 
The dephasing rate we used in calculation T2/2n — 0.05GHz, 
and the intra-well relaxation rate rio/27r = 0.6GHz. 



where W« is the LZ transition rate introduced in Eq.©, 
and Tij is the relaxation rate from \i) to \ It is noticed 
that some relaxation rates (e.g., Toi) that have little ef- 
fects on the results have been neglected for simplicity , 
and these rates would be considered in part C. Moreover, 
we did not consider the population of |3) for two reasons: 
(i) The driving-field amplitude is not large enough to 
reach anti-crossing A03, resulting in no LZ transitions 
between |0) and |3). (ii) In general, for the supercon- 
ducting flux qubits, the intra-well relaxation rate Tio is 
much larger than the LZ rate W13. Therefore, the pop- 
ulation of 1 1) mostly relaxed to |0) rather than made a 
transition to |3). In part C these approximations would 
be discussed in detail and we found they are valid for the 
reported experimental data. 

In the stationary case, t>o = Pi = p 2 = 0. The qubit 
population in the left well can be easily solved from 
Eq.©: 



P2 



W 02 (W 12 +F 10 ) 



w 12 (2W 12 + w 02 + r 20 ) + r 10 (2W 02 + w 12 



r 20 ) 

(8) 

Fig. 2(c) is the contour plot according to Eq.(T4]) and 
Eq.© with the parameters from experiments 2 ^. The 
main features of the theoretical results agree with the 
experimental data very well. Since the pattern is sym- 
metric about detuning, we only plot half of the pattern. 
Qualitatively, using Eq.©, we split the space into three 



FIG. 3: (color online) (a) LZ transition rate W03 versus flux 
detuning and driving amplitude. (b)LZ transition rate W12 
versus flux detuning and driving amplitude, (c) Qubit pop- 
ulation in the left well obtained from Eg, (11011 . Features of 
population inversion and checkerboard pattern are notable. 
The parameters we used are the same with those in Fig. 2. 



regimes: 



P2 



W 02 



W02 is off 

■>u I-.,, w 02 is on;Wi 2 is off 

$£->0 W 02 is on; W 12 is on. 



When the amplitude of rf field is small, Wq 2 is off. 
The channel that connects states |0) and |2) is cut off, 
resulting in no population transfer to the left well. The 
minimum amplitude required to switch Wq 2 on is propor- 
tional to the detuning. Therefore, we observed an edge 
with positive slope at left. With the amplitude increased, 
we may enter a regime where Wq 2 is on and Wi 2 is off. 
The qubit serves as a two level system with the chan- 
nel to higher energy levels cut off. Since the minimum 
amplitude to switch on the W12 linearly decreases with 
the detuning, we saw a shadow with negative slope on 
the right side. There is another regime where both Wq 2 
and W\ 2 are on. Because W12 is much larger than W02, 
the population of |1) is more likely to be excited to |2), 
and then relaxes back to |0), thus suppressing the net 
population transfer to the left well (see the red path in 
Fig.l). In addition, the population on the peaks of the 
interference patterns should be less than 0.5 because of 
the spontaneous relaxation to the ground state. This was 
confirmed in the experiments. 



B. THE SECOND DIAMOND 

The 'second-diamond' interference patterns correspond 
to the transition processes shown in Fig.l (green path). 
Using a similar method with that of part A and neglect- 
ing weak transitions, we can write the rate equation for 
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FIG. 4: (color online) Calculated qubit population versus flux 
detuning and driving amplitude. The first diamond corre- 
sponds to the red path and the second diamond corresponds 
to the green path in Fig.l. The parameters we used in calcu- 
lation are: the inter-well relaxation rate T2o/2tt = 0.05MHz, 
the avoided crossing Ai3=0.5GHz, and the temperature T 
used in Fq2 = F2o exp(—Eo2/k b T) is 20mK. Other parame- 
ters are the same with those in Fig. 2. 



the population evolution as: 



Po = -PoW 03 +PiT 1Q +p 2 r 2 o 
Pi = -pi(Wn + Tio) + p 2 Wi 2l 

pi =piW 12 -P2{Wi2 + r 20 ) +P3r 3 2, 

PO +P1+P2 +P3 = 1, 



03: 



(9) 



where Wij and T^ have the same definitions as those in 
part A. In the stationary case po = Pi = P2 = Po — 
0, Eq.© can be analytically solved. Considering Tio, 
^> Wi2,Wo3> which is the case in experiments, the 
population in the left well can be obtained in a very sim- 
ple form: 



PL =P2+ P3 



in 



Wo3 + w 12 + r 2 



(10) 



Fig. 3(c) is the contour plot using Eq. (fl0|) with exper- 
imental parameters. The results, which show notable 
features of both checkerboard patterns and population 
inversion, agree with those of the experiments very well. 
The population inversion is a straightforward result of 
Eq. (10). When W 03 is on and W%2 is off, since 
Wq3 > T 2 o, the system forms a A- type three- level struc- 
ture, pl is larger than 0.5 and even close to 1. Actually, 
one can use this property to generate controllable pop- 
ulation inversion to produce lasers. However, once W12 
is on, the population of \2) is pumped to |1) and relaxes 
to |0) rapidly. Because W12 > W03, pl < 0.5, thus de- 
stroying the population inversion(see the green path in 
Fig.l). 

C. COMBINE THE TWO DIAMONDS 



FIG. 5: (color online) (a) Qubit population obtained from 
Eg . (111 [I with driving frequency ui/2n — 1.2GHz. Moire-like 
patterns reveal for resonances n > 12. The parameters we 
used are from experiments»*ii, with A02 = 0.004GHz and 
A12 is fitted to be 0.09GHz. The locations of the level 
crossings A02 and A12 are OrmFo and 13.1m"l>o respectively. 
The diabatic energy-level slope |mo|(|m2|)=1.01GHz/m<I>o, 
and |mi|(|m3|)=0.91GHz/m<I , o. The dephasing rate we used 
in calculation r2/27r = 0.05GHz, the inter- well relaxation 
rate r2o/27r = 0.05 MHz, and the intra-well relaxation rate 
Fio/27r = 0.6GHz. (b) The dephasing rate is increased to 
r2/2-7r = 0.2GHz, leading to a diamond-like interference pat- 
tern, in which the area of missing parts increases. 



Having addressed the two diamonds separately, now we 
can make a more integral and quantitative study of the 
interference patterns. Here the driving-field amplitude 
is increased from to a large value which has reached 
Arj3 and A 12 crossings. In this situation, all the transi- 
tion processes in Fig.l should be considered, and the rate 
equation can be written as: 

Po = -Po{W 2 + W03 + T 02 ) +piT la +p 2 (W 2 o + r 20 ) 
+ P3W03, 

Pi = -pi(W 12 + W13 + T w ) + p 2 W 12 + P3VF13, 

t>2 = po{w m + r 02 ) + P1W12 -P2(w Q2 + W12 + r 20 ) 

+ P3^32, 

P0+Pl+P2+P3 = l, (11) 

which can also be solved analytically. Since the analyt- 
ical solution is too complex to extract a physical picture, 
we did not write it out explicitly here. Nevertheless, we 
solved Eq. (jlip numerically and obtained the qubit pop- 
ulation as functions of the amplitude and flux detuning as 
shown in Fig. 4. The agreement between our results and 
those of the experiments is remarkable. The right edges 
of the second diamond are not shown in Fig. 4 because 
it requires the presence of higher excited states, which 
could be considered easily with the method we discussed 
here. 

In addition, there is no explicit signature of coher- 
ent traversal through avoided crossing A13 observed in 
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FIG. 6: (color online) Qubit population obtained from Eq, (|ll[) with driving frequency u>/2n = 0.16GHz, 0.879GHz and 
0. 8886GHz from left to right, and dephasing rate r2/2-7r=0.05GHz, 0.2GHz and 0.4GHz from top to bottom, respectively. The 
parameters of the superconducting flux qubit are the same with those in Fig. 2. 



experiments—. In order to clarify this point, we changed 
the value of A13 from to 0.8GHz. The interference pat- 
terns obtained from Eq.(ll) exhibit almost no change, 
indicating that the LZ transitions at A13 has negligible 
contribution to the whole process. Moreover, we have 
compared the results of part C to those of part A and 
B respectively. They accorded very well, indicating that 
our approximation in part A and B is appropriate. 

III. EFFECT OF DRIVING-FREQUENCY 
AND DEPHASING RATE 

As discussed in ref [12], subsequent LZ-tunncling 
events can interfere only when phase coherence is pre- 
served, requiring that the time interval of subsequent 
LZ tunneling events t p must be less than de-coherence 
time T co h of the system (t p < r co h ). Since r v oc lo^ 1 , if 
the qubit 's decoherence time T co h is short due to strong 
coupling to environment, a large driving frequency is re- 
quired to observe LZ interference. Therefore, it is of great 
meaning to discuss the behavior of the qubit driven by 
high-frequency fields. We emphasize that the high fre- 
quency we discussed here is still in the regime dominated 
by LZ transitions, which means that the analysis in Sec. II 
is still valid. 

Another advantage of using high frequency is that one 
can resolve more features in the interference patterns. At 



low frequency, the microwave can match the energy level 
spacing approximately at all flux detuning, correspond- 
ing to the almost continuous band of Wij in Fig. 2 and 
Fig. 3. Actually, these continuous bands are composed 
of resonant peaks (with single and multiple photon ef- 
fect) that are overlapped with each other. The distance 
between these peaks is the driving frequency. At high 
frequency, the distance of these overlapped resonances 
increased and the continuous band in Fig 2. and Fig. 
3 become a series of discrete peaks. This was demon- 
strated by recent experiments. By using higher driving 
frequency, clear LZ interference pattern is observed and 
the results show some new characters(e.g., partly missing 
fringes) that are different from the 'diamond' interference 
patterns^. We found that L2 ±2 u/2ir is the critical sit- 
uation that the interference patterns transform between 
moire-like and diamond- like types. To explain the moire- 
like pattern^ observed for resonances n > 12, where n 
represents the n-photon resonances, higher excited states 
need to be considered. By inserting the experimental 
parameter s 9 ! 11 into Eq.(ll), the qubit population was ob- 
tained (see Fig. 5(a)). The calculated results, with some 
fringes partly missing and some fringes intact, agree with 
the experiments very well^ These features can be eas- 
ily understood using Eq. (9). The missing part of the 
fringes at n > 12 corresponds where both W02 and W12 



6 



are on. The intact fringes (e.g., n=14) correspond where 
W 02 is on while Wi 2 is off. 

For higher dephasing rate, the width of the resonance 
peaks for all Wij will increase. As a result, the area of 
the region where both W02 and Wn are on also expands, 
resulting in an increasing area of missing fringes(see Fig. 
5(b)). Similar to that of the part A, the missing fringes 
form a shadow with negative slope. 

We investigated the effect of driving frequency and 
dephasing rate on the LZ interference. Shown in fig. 6 
are some examples. With dephasing rate increasing, the 
individual resonances are no longer distinguishable and 
merged into a continuous band (Fig. 6 (d-f)). However, 
by increasing the driving frequency to w/2tt > T2, we re- 
build the interference again (Fig.6(d,g)), as discussed in 
ref[13]. It is interesting that there are some new features 
emerging during these processes. 

As we discussed in Sec. II, whenever W03 and W12 are 
on simultaneously, the population will be pumped back 
and a missing part will be generated on the fringes (Fig. 
3). For small dephasing rates, the resonant peaks are 
very sharp. The resonant peaks are difficult to hit the 
spots where both W03 and W\i are on. We thereby saw 
high quality moire-like interference patterns. However, 
for large dephasing rates, the resonant peaks are very 
wide along the axis of flux detuning. At many places 
W03 and W\i manipulate the resonant peaks simulta- 
neously and we observed complicate patterns that are 
mainly diamond-like. Similar anomalous pattern was re- 
ported in recent experiments^. 

IV. CONCLUSION 

We have analyzed the interference patterns observed 
in recent experiments which demonstrated an innovative 



approach to make spectroscopic measurement of a quan- 
tum system. The analytical results obtained from our 
model are agreed with the experiments very well, and 
the method can be extended to investigate the response 
of higher energy levels under a large-amplitude driving 
field. Moreover, we investigated the effect of driving fre- 
quency and dephasing rate on the LZ interference. The 
fast dephasing will destroy the interference. However, it 
is possible to overcome the dephasing by increasing rf fre- 
quency. We have not only explained the moire-like pat- 
tern observed in a recent experiment, but also discussed 
some interesting anomalous patterns that could be ex- 
pected in experiments. Our results can be used to un- 
derstand the LZ interference in multilevel system under 
different driving frequencies and decoherence rates. The 
model and method used in this article can be extended 
to other systems with multiple energy levels structure. 

In addition, we noticed a recent work-22 in which 
microscopic two-level systems (TLS) were observed in 
flux qubits by using one- and two-photon spectroscopy 
method. We proposed that it is very convenient to detect 
such TLSs with the amplitude spectroscopy method 
demonstrated by Berns et al^ and the analytical 
methods in this article will be still valid in such case. 
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